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Abstract 

Motivated by the importance and universal character of phase singular- 
ities which are clarified recently, we study the local structure of equi-phase 
loci near the dislocation locus of complex valued planar and spatial waves, 
from the viewpoint of singularity theory of differentiable mappings, initi- 
ated by H. Whitney and R. Thom. The classification of phase-singularities 
are reduced to the classification of planar curves by radial transformations 
due to the theory of A. du Plessis, T. Gaffney and L. Wilson. Then fold 
singularities are classified into hyperbolic and elliptic singularities. We 
show that the elliptic singularities are never realized by any Helmholtz 
waves, while the hyperbolic singularities are realized in fact. Moreover, 
the classification and realizability of Whitney's cusp, as well as its bifur- 
cation problem are considered in order to explain the three points bifurca- 
tion of phase singularities. In this paper, we treat the dislocation of linear 
waves mainly, developing the basic and universal method, the method of 
jets and transversality, which is applicable also to non-linear waves. 

1 Introduction. 

A complex scalar wave has the locus, the dislocation locus, where its phase is 
not denned. The local structure of equi-phase loci near the dislocation locus is 
called a phase- singularity |21| . The phase singularities are called optical vortices 
in optics and are very basic and important objects in any science related to 
waves and quanta. In this paper we give the exhaustive classification of phase 
singularities of complex scalar waves of low codimension. 

In [2nj|2J|, J- F. Nye constructed extensively complex scalar global planar 
waves satisfying the Helmholtz equation, with detailed analysis of those exam- 
ples. Also he gave, by his examples, an explanation of an experimental bifur- 
cation process of phase singularities: one degenerate singular point bifurcates 
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to three singular points and then another singular point annihilates with one of 
the three. He intends to explore the phase singularities from an analogy with 
the catastrophe theory 

In this paper, we understand phase singularities clearly using the singularity 
theory of differentiable mappings [3] JU| E3 HI ■ 

The planar complex scalar wave can be regarded, from the general point of 
view, simply as a differentiable mapping from the plane to the plane of complex 
numbers. Then, by a theorem of H. Whitney the generic singularities of the 
wave, as a differentiable mapping, are just the fold singularities and the cusp 
singularities. The singular values form on the plane of complex numbers an im- 
mersed curve, the discriminant, with several number of cusps. Then generically 
the discriminant does not hit the zero, so that the zero is a regular value; generic 
phase singularities are regular, namely, locally diffeomorphic to the standard ra- 
dial lines emitted from the origin. However, for a generic time-depending wave, 
the curve of singular values moves and momentarily may hit the zero. Thus, 
generically momentary wave can have degenerate phase singularities described 
by the fold singularities. Moreover, for a generic two parameter family of plane 
wave, the cusp singularity occurs as more complicated phase singularities. 

The above simplified story must be examined twofold: First, in Whitney's 
theorem, the singularities are classified by means of arbitrary local diffeomor- 
phisms of the source plane and the target plane. However, for the classification 
of phase singularities, we concern with the equi-phase lines and thus need to 
consider finer classification using only diffeomorphisms which preserve the ra- 
dial lines on the target plane. Second, because waves must obey several natural 
conditions given by, say, the Hclmholtz equations and the wave equations, more 
than just the differentiability, we must consider the realizability of singularities 
and determine generic singularities among waves satisfying those conditions. 

We clarify the equivalence relations for phase singularities, and thus classify 
all phase singularities of low codimension, and discuss the realizability by the 
Helmholtz waves of phase singularities. Further, we propose the new explanation 
for the experimental bifurcation process treated in |2U | |21 ) . 

In the next section, we formulate our equivalence relation providing the base 
of our classification. A natural and refined classification by radial transforma- 
tions is established on phase singularities for planar and spatial complex scalar 
waves. Then, we give the exact classification of generic complex planar waves 
without conditions motivated from physics. 

In §3, the realization of singularities by the Helmholtz waves is examined 
by concrete examples, which have a different character with Nye's examples in 

In §4, we treat phase singularities of spatial complex scalar waves and con- 
sider their realizability. 

In ^S] the classification problem of planar waves is reduced to that of planar 
curves under diffeomorphisms preserving radial lines. 

We introduce in ^jthc notion of Helmholtz jet spaces and transversality to 
discuss genericity of singularities for Hclmholtz waves. 
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In 33 as an application of the method developed in this paper, we discuss the 
bifurcation problem of phase singularities of solutions to non-linear Schrodinger 
equations. 

In this paper, we consider local classification problem of phase singularities. 
For the global topology of dislocation locus, see 

For other applications of the singularity theory to solutions of partial differ- 
ential equations, see ms ^ ance - 

2 Phase singularities for planar complex scalar 
waves. 

We denote by C the plane of complex numbers and write a complex number 
as u + iw = re 10 , u, w being the real part and the imaginary part respectively, 
while r, 9 the modulus (or the amplitude) and the argument (or the phase) 
respectively. 

Let us consider a complex scalar wave 



on the (x,y)-plane depending on the time (or any other one-parameter). First 
we regard 'J/ as just a time-depending complex valued function. We assume 
u(x,y,t), v{x,y 1 t) are differentiable (i.e. C°°) functions. 

If ^>(x,y,t) =/= at a point (x,y) and at a moment t, then we can write 
^(x, y, t) = r(x, y, t)e l9 ( x,v,t ^ uniquely with r(x, y, t) > and 8(x, y, t) mod. 2tt. 
Then, we are concerned with the wave dislocation locus at a moment t = to 

{(x, y) | 1>(x, y, t ) = 0} = {(x, y) | u(x, y, t ) = 0, w(x, y, t ) = 0} 

and the equi-phasc curves {(x,y) \ 9(x,y,to) — const.} outside of the wave 
dislocation locus. 

Then, in the framework of singularity theory of differentiable mappings, we 
introduce the notion of radial transformations and give the exact classification 
results of singularities relatively to the radial transformations. 

A radial transformation on C near is a diffeomorphism, an invcrtible dif- 
ferentiable transformation, t(u,w) = (U, W), r: (C,0) — > (C,0) which sends 
any radial line {8 = const.} to a radial line. In fact, a diffeomorphism t(u, w) — 
(U, W) is a radial transformation if and only if there exists a positive function 
p(u, w) and real numbers a, 6, c, d with ad — be =/= such that 



For the classification, we define the equivalence relation on phase singu- 
larities: Two functions ^(x, y,t) and ^{x,y,i) — u'{x,y,t) + iw'{x, y,t) are 
radially equivalent at points and moments (xo,yo,to) and (x' ,y ,t' ) respec- 
tively, if there exist a local diffeomorphism a(x,y) = (X(x,y),Y(x,y)) on the 



* = 9{x,y,t) 



u(x,y,t) +iw(x,y,t) 



U = p(u, w)(au + bw), 



W = p(u, w)(cu + dw). 
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plane with X(xo,yo) — x' ,Y(xo,yo) = y' and a local radial transformation 
t{u, w) — (U(u, w), W(u, w)) near the origin on C such that 

u(X(x,y),Y(x,y),t ) = U(u'(x,y,t' Q ),w' (x,y,t' Q )), 
w(X(x,y),Y(x,y),t Q ) = W(u' (x,y,t' ),w' (x,y,t' )), 

namely that $(a(x,y),t ) = T($(x,y,t' )). 

Theorem 2.1. For a generic complex valued function ^(x,y,t), the phase sin- 
gularity at any point and any moment (xo,yo,to) is equivalent under radial 
transformations to the regular singularity 

R: if)(x,y) = x + iy, 

the hyperbolic singularity 

H: ip(x,y) = x 2 - y 2 + iy, 
or to the elliptic singularity 

E: ip(x,y) = x 2 + y 2 + iy, 
at the origin (x,y) = (0,0). (see FigureQJ ) 

Each phase singularity of the classification in Theorem 12. II is determined by 
its two jet actually. 




R H E 



Figure 1: phase singularities 

Remark 2.2. Besides the wave dislocation, we can classify generic critical 
points of phase functions defined outside of the dislocation locus: The generic 
critical points are the non-degenerate maximal( minimal) points, the saddle points 
and the cuspidal points. The last one bifurcates to one maximal(minimal) point 
and one saddle point . 

Both hyperbolic and elliptic singularities are equivalent to the fold singularity 
if) : (x,y) i-> (u,w) = (x 2 ,y) 
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under arbitrary diffeomorphisms not necessarily radial transformation, namely, 
under the right-left equivalence. 

For a momentary complex wave tp(x,y) = u(x,y) + iw(x,y) on the plane, 
the locus in C of complex values ip(xo,yo) for (xo,yo) with 



di 



du 



dot 



dx dy 
dw dw 

dx dy 



(£o,2/o) = 0, 



is called the discriminant of the complex wave ip- 

Note that the above normal form of fold singularity is never generic as the 
phase singularity. In fact, the discriminant of ip in that case is the w-axis 
in C which has the infinite tangency (actually coincides) with the radial lines 
{9 = tt/2} and {9 = 37r/2}. Generically the discriminant must be tangent to the 
radial lines in non-degenerate manner, namely, in the second order tangency. 
Then there are two possibility of non-degenerate tangency of the discriminant 
of fold singularities at £ C; the image of ip (the value set of ip) is concave 
or convex. These correspond, respectively, to the hyperbolic singularity and 
elliptic singularity. 

Moreover, the generic bifurcations on t of the hyperbolic singularities and 
the elliptic singularities are given by 



H t : ®(x,y,t) 
E t : *(x,y,t) 



x 2 ~y 2 +t + iy, (teR), 
x 2 + y 2 + t + iy, (teR). 



(see Figure El) 







t<0 



t=0 



t>0 



Figure 2: The bifurcations of the hyperbolic phase singularity (top) and the 
elliptic phase singularity (bottom). 
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The picture of the bifurcation of the hyperbolic singularity can be seen in 
Fig. 6 of 0. 

Remark 2.3. The classification of phase singularities is closely related to the 
classification of plane curves under diffeomorphisms preserving a given singular 
foliation on the plane. Then, one of the most delicate cases is the case when 
the foliation is formed by radial lines, which is given by Euler vector held X = 
d d 

x— — h y— (cf. That is the case we are treating in this paper (see §5). 

ox ay 

The discriminant of the fold singularity is a regular curve. Then degenerate 
phase singularities are classified as follows: 

Proposition 2.4. The phase singularities arising from fold singularities are 
classified into 

ip m (x,y) = x 2 ±y m + iy, to = 2, 3, 4,..., 

under radial transformations (and diffeomorphisms on the target), provided the 
discriminant curve has a contact with the tangent line at the origin in a finite 
multiplicity. 

Remark 2.5. The generic bifurcation of il> m (x,y) is described by the family 

x 2 ± y m + t m -iy m ~ l + t m ^ 2 y m - 2 + ■■■ + t 2 y 2 + t + iy, 
with (to — l)-parameters to 1*2, • ■ ■ , t m _i. 

A momentary complex wave ip : (R 2 , (xq, yo)) — > (C, 0) is called a Whitney's 
cusp or simply a cusp if it is right-left equivalent (under local diffeomorphisms on 
R 2 and C which are not necessarily radial) to the mapping ip(x, y) = x 3 +xy+iy. 
We are interested in this type of phase singularity because there occurs a three 
points bifurcation by just a translation ip a (x, y) — x 3 +xy+i(y+a), (a £ R). (For 
the classification of more degenerate singularities under the right-left equivalence 
relations, see [22] EH])- 

The Whitney's cusp appears generically in two parameter families of planar 
complex valued functions. 

Then, we have 

Proposition 2.6. (The radial classification of Whitney's cusps): Any Whit- 
ney's cusp is equivalent under radial transformations to the standard function 
ip(x, y) = x 3 + xy + iy. 

The typical bifurcation of the phase singularities for a Whitney's cusp is 
described by *p a ,b(x, y) = x 3 + xy + b + i(y + a), (a, b S R). (see FigureOJ 

Remark 2.7. The bifurcation problem of phase singularities arising from Whit- 
ney cusps is related to web geometry (pQ, ^0]). In fact, generic two parameter 
families of Whitney cusps define 3- webs on the plane, and their classification by 
radial transformations provides functional moduli (R,cmark l5.4|l . 
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Figure 3: The two parameter bifurcation of the cusp phase singularity. 

3 Phase singularities of the Helmholtz waves. 

Now, we ask the physical reality of the classification; the instantaneous appear- 
ance of singularities for wave functions satisfying the wave equation and the 
Helmholtz equation. Namely, we assume the wave ^(x,y,t) satisfies the wave 
equation 

d 2 <S> 



dt 2 



c 2 V 2 *, 



d 2 d 2 

for a positive real number c, where V 2 = -7—^ + — ■= is the Laplacian. Moreover 

ox z ay 1 

we assume that ^(x,y,t) satisfies the Helmholtz equation 

V^ + Zc 2 * = 0, 

for a positive real number k as a very natural physical assumption for monochro- 
matic waves. We call a function which satisfies the Helmholtz equation the 
Helmholtz function. Further, we call the Helmholtz function with a parameter 
t which satisfies the wave equation as well the Helmholtz wave. Note that if we 
have a solution ^ for c = 1, k = 1, then, by setting fy(x, y, t) = ^{kx,ky, £t), 
we have a solution 'J for general c and k. 



7 



By solving the Cauchy problem properly, we obtain the following. 

Proposition 3.1. The complex valued function 

ip(x, y) = cosy — cos 2: + isiny 

satisfies Helmholtz equation ip xx + tp yy + ip — and has the hyperbolic singularity 
at the origin. Moreover the hyperbolic singularity with its generic bifurcation is 
realized by a Helmholtz wave 

^(x, y, t) = (cosy — cosx + isiny) cost + cos y sint, 

(for k = 1, c = 1). 

Wc sec it is radially equivalent to the normal form simply by observing its 
Taylor expansion. 

To the contrary, we observe: 

Proposition 3.2. Elliptic singularities are not realized as a function satisfying 
the Helmholtz equation. 

Proof: Suppose a function ^(x^y) is radially equivalent to the elliptic singular- 
ity. Then, the image of ip is convex at the origin where the tangent line supports. 
Suppose the function ip(x, V) satisfies the Helmholtz equation with k = 1. From 
the equation tp xx + ip yy + tp — 0, we see the Hessian of tjj is traceless at the 
dislocation locus {ip = 0}, so are the real part Re(HessV') and the imaginary 
part Im(HessV'). Thus, for any real numbers A, /1, 

ARc(HessV') + /xlm(HessV') 

is never a definite matrix. However, the linear projection along the tangent line 
to the image of ip must be definite. This leads to a contradiction. □ 

In fact, as for the generic classification of phase singularities for one-parameter 
families of complex valued functions satisfying Helmholtz equation, we have: 

Theorem 3.3. The generic phase singularities of planar Helmholtz functions 
are regular singularities and hyperbolic singularities. 

Theorem 3.4. The generic phase singularities of planar Helmholtz waves are 
regular singularities and hyperbolic singularities. 

Theorem 13.31 and Theorem 13.41 are proved in Sectional 

For the cusp singularities, we have: 

Proposition 3.5. A Whitney's cusp is realized as a Helmholtz wave. In fact, 

ip(x,y) = x 3 cos y + (x — 3xy) sin y + i sin y 

is a Whitney's cusp satisfying Helmholtz equation (k — 1): ip xx + ip yy +-0 = 0. 
Moreover, 

*fr(x, y, t) = (x 3 cos y + (x — 3xy) siny + isiny) cost + i cosy sin t 

gives a deformation of ip by a Helmholtz wave (k = l,c = 1) describing a three 
point bifurcation of the phase singularity. 
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Remark 3.6. By a similar construction to Proposition 13.51 we have another 
realization 

i(>(x, y) = x 2 cos y — y sin y + i sin y 
of hyperbolic singularities. 

Remark 3.7. Apart from the classification problem of phase singularities, we 
can show that generic Helmholtz function ip : R 2 — > C is, locally at any point in 
R 2 , right-left equivalent to a regular point, to a fold point or to a cusp point. 

4 The radial classification of phase singularities 
for spatial complex scalar waves. 

We study, in this section, the phase singularities of spatial waves \& = ^(x, y, z, t) : 
R 3 x R — > C. 

The generic singularities of differentiable mappings R 3 — > C consist of the 
definite fold singularities, the indefinite fold singularities and the cusp singular- 
ities. The normal forms of them are given by 

the definite fold singularity : ip(x, y, z) — x 2 + y 2 + iz, 
the indefinite fold singularity : ip(x,y,z) — x 2 —y 2 + iz, 

the cusp singularity : ip(x,y,z) — x 3 +xy + z 2 + iy, 

under the left-right equivalence |14|. For a generic complex valued function 
^{x, y, z, t), only fold singularities may appear as a phase singularity. Moreover 
the discriminant curve has non-degenerate tangency with the tangent line at 
the origin. Thus we have 

Theorem 4.1. For a generic spatial complex valued function ^(x,y,z,t), the 
phase singularity at any point and any moment (xq, j/o, Zq, to) is equivalent, un- 
der the radial transformation, to the regular singularity 

R: ip(x, y, z) = x + iy, 

to the definite hyperbolic singularity 

DH: ip(x, y, z) = x 2 + y 2 — z 2 + iz, 

to the definite elliptic singularity 

DE: ip(x, y, z) = x 2 + y 2 + z 2 + iz, 

or to the indefinite singularity 

I: ip(x, y, z) — x 2 — y 2 — z 2 + iz, 

at the origin (x,y,z) = (0,0,0). (See FigureQJ) 
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DE 
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Figure 4: Phase singularities for spatial waves 



Similarly to the case of planar complex scalar waves, the generic bifurcations 
on t of the definite hyperbolic singularities, the definite elliptic singularities, and 
the indefinite singularities are given by 



DH t 
DE t 
It 



y>{x, y,z,t) = x 2 + y 2 -z 2 +t + iz, 
V(x, y,z,t) = x 2 +y 2 + z 2 + t + iz, 
^{x, y, z, t) — x 2 — y 2 — z 2 + t + iz. 



(see Figure 0) 



DHt 





t<o 



DEt 








t 


=0 




















t<0 



t=0 



It 





t<0 t=0 t>0 

Figure 5: The bifurcations of phase singularities of spatial scalar waves. 

For cusp singularities, we have 

Proposition 4.2. The phase singularities which come from cusp singularities 
are all radially equivalent to 



i/j(x, y, z) — x 3 + xy + z 2 + iy. 



10 



Proposition 14. 21 is induced from Lemma \b. 31 

A complex valued function ^(x, y, z,t) is called a Helmholtz wave if it sat- 
isfies the wave equation ^> t t = c 2 ($„ + ^ yy + $ zz ) for a positive real number 
c and the Helmholtz equation ty xx + ^ yy + ^ zz + k 2 ^ = for a positive real 
number k. 

As for the realizability of the spatial waves as Helmholtz waves, we have: 

Proposition 4.3. The definite hyperbolic singularity and the indefinite singu- 
larity together with their generic bifurcations are realized by Helmholtz waves 
(for k = 1, c = 1): 

DH t : y, z, t) = (— cos x — cos y + 2 cos z + i sin z) cost + cos z sint. 

It : ^(x, y, z, t) = (— 2 cosa; + cosy + cos z + isinz) cost + coszsint. 

Moreover, in a similar way to Proposition 13. 21 we have 

Proposition 4.4. Definite elliptic singularities are not realized as Helmholtz 
waves. 

The cusp singularity is realized as a Helmholtz function: 

Proposition 4.5. The complex valued function 

ip(x,y, z) — x 3 cos y + (x — 3xy) sin y — cos y + cos z + i sin y 

satisfies Helmholtz equation ip xx + ip yy + ip zx +-0 = and is radially equivalent 
to the cusp singularities. 

5 Radial classification of planar curves. 

The classification problem of complex waves under radial transformations is 
reduced to the classification problem of planar curves under radial transforma- 
tions, by means of du Plessis, Gaffney, Wilson's theory [TI][T3]. The theory 
reduces the classification to that of discriminants with an exception. The ex- 
ceptional cases are hyperbolic and elliptic singularities. Although they have the 
same discriminants, they are not radially equivalent. It depends on whether the 
image is convex or concave. The following results on curves have been applied 
to the classification of discriminants, and then, phase singularities. 

Lemma 5.1. (The radial classification of regular curves): A regular curve 
through the origin on C is transformed by radial transformations to the curve 
u = w m for some integer m > 2 or u — g(w) for some function with null 
derivatives g (i) (0) = 0,i = 0, 1, 2, 3, . . . . 

Proof: Using a linear transformation, we may suppose the regular curve is given 
by u = g(w) for a function g(w). Suppose ordg = at w = and <r m )(0) > 0. 
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Then we can write u = (a(w)w) m for a function a(w) with a(0) > 0. Set 
p{w) = a(w) m ~ 1 , and define the radial transformation U — p(w)u, W = p(w)w. 
Then U = p{w){a(w)W/p(w)} m = W m . □ 
In Section [21 we study deformations of phase singularities. For them we 
observe the following. 

Remark 5.2. Let (u(t, X),w(t, A)) be a deformation of the curve (u(t, 0), w(t, 0)) = 
(t m ,t): 

u(t,X) = ao(A) + ai(X)t H \-a m (X)t m -\ , 

w(t,X) = /3 (A)+/3i(A)t + --- , 

with ao(0) = • • • = a m -i(0) = /3 (0) = 0, a m (0) = /3i(0) = 1. Then by a family 
of radial transformations, the family is transformed to 

u(t, X) = a' (X) + a' 2 {X)t 2 + ■■■+ ^(A)*™- 1 + t m , 
w(t,X)=l3' (X)+t, 

for some functions a' (X) , a' 2 (X) , . . . , a' m _ 1 (X) , 0' O (X) . The latter curve is ex- 
pressed as 

u = t (A) + ii(A)ic + • ■ ■ + t m ^i(X)w m ~ 1 + w m . 
By a family of linear transformations (u, w) i— > (u — ti(A)u>, w), it is reduced to 

7i = to(A) + t 2 (X)w 2 + ■■■+ t m ^ 1 (X)w m - 1 + w m . 



In general, any parametrized curve (u(t),w(t)) through the origin in C is 
equivalent by radial transformations and re-parametrizations to 

u(t) =t m + 0(t m+1 ), w(t)=t n , 

for some integers m, n with m > n. We have n = 1 for regular curves. If n > 2, 
then we call the curve an (n,m)-cusp. A (2,3)-cusp is called a simple cusp 
briefly a cusp. 

Lemma 5.3. (The radial classification of simple cusps): Any simple cusp is 
equivalent by radial transformations and re-parametrizations to 

u(t)=t 3 , w(t)=t 2 . 

Thus, any two simple cusps are radially equivalent to each other. 

Lemma 15.31 implies Proposition 14.21 since the discriminant ofip(x,y,z) — 
x 3 + xy + z 2 + ly is a simple cusp. 

Proof of Lemma \5.3\ : Let u(t) = t 3 + 0(t 4 ), w(t) — t 2 be a simple cusp. By a 
linear transformation on the (u, w)-plane and a re-parametrization of t, we may 
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suppose the curve is given by u(t) = t 3 , w(t) = t 2 + 0(t 4 ). Set w(t) = t 2 a(t) for 
a smooth function a(t). Then a(0) = l,a'(0) = 0. Then there exists a smooth 
function p(x,y) such that a(t) — p(t 2 ,t 3 ) by the preparation theorem (HU). 
Then the curve is radially equivalent to the curve 

«<*> = ^Ttaja* 8 ' = ^T^' 

which is radially equivalent to u(t) = t 2 ,w(t) = t 3 . □ 

Remark 5.4. Let C a ,b(t) — {u(t, a,b),w(t,a,b)) be a generic two parameter 
family of simple cusps. For each (a, b), we draw tangent lines to the simple 
cusp C a .f, from the origin. Then there exists a non-void open subset U such 
that for (a, b) E J7, there are exactly three tangent rays. By the assignment 
of the corresponding tangent points, we have three functions Ai, A2, A3 on 
U ; C 0) &(Ai), C a ,i>(A2), C a ,i(A3) are tangent points. Thus we have a triple of 
foliations: 

Ai = const., A2 = const., A3 = const., 

that is, a 3- web on U. Moreover, radially equivalent families of simple cusps have 
isomorphic 3-webs. It is known that the classification of 3-webs has function 
moduli in general (|1(J)'). 

6 Helmholtz jet space and transversality. 

We introduce the notion of Helmholtz jet spaces and show the transversality 
theorem in a Helmholtz jet space, as one of the main ideas to show the results 
in this paper. Note that, in JS], analogous jet spaces are considered for other 
kinds of Monge- Ampere equations. 

Consider the Taylor expansion around (x, y) = (xo,yo) of a complex valued 
function ip on the (a;, y)-plane: 

tl>(x,y) = a+bX+cY+lx 2 + fXY+ 9 -Y 2 + ^X 3 + ^X 2 Y+ t -XY 2 + ^Y 3 + --- . 

I I 2 2 D 

Here we set X = x — Xq, Y = y — yo, and a, 6, c, . . . are complex numbers. 

Suppose ip is a Helmholtz function for k = 1, that is, "0 satisfies Helmholtz 
equation ip xx + ijj vy + ip = 0. Then we have 

e + g + a = 0, h + £ + b = 0, k + m + c = 0. 

Therefore, we have 

ip(x,y) = a + bX + cY + ^X 2 + fXY -^(a + e)Y 2 

+ ^X 3 + \x 2 Y \{b + h)XY 2 -l(c+ k)Y 3 + ■■■. 

The Taylor expansion of a function ip up to order r around a point (xq, yo) 
of R 2 is called the r-jet of ip at (xo,yo) and denoted by j r ip(xo,yo)- Denote 



13 



by J r (R 2 , C) the space of r-jets of complex valued functions on R 2 . In it, we 



denote by J^ clm (R , C) the set of r-jets of planar Helmholtz functions for k = 1: 
,2 ,C) = {fi>{x ,y ) | ip xx + ipyy + ip = around (x ,y )}. 



We call it the Helmholtz r-jet space. For example, JH olm (R 2 ,C) is identified 
with R 16 = R 2 x C x C 6 with coordinates xq, yo, a = a± + i<i2; b = b\ + i& 2 , 
c = ci + ic 2 , e = ei + ie 2 , f = fi + i/2, h = h% + ih 2 , and fc = fci + ifc 2 . 

In general, the Taylor expansion of a Helmholtz function ip(x,y) defined 
around (xo,yo) is determined by ip(x,y ) and ip y (x,yo). Moreover, for any 
given complex valued analytic functions ipo(x) and ipi(x) defined around Xo, 
there exists uniquely a complex valued function ip(x, y) defined around {xq, yo) 
satisfying the Helmholtz equation, ip(x, yo) — ipo(x) and ip y (x, Vo) = ipi{x). The 
r-jet of ip at (xq, yo) is determined by the r-jet of ipo(x) at xq and (r — l)-jet of 
ipi{x) at xq. Thus, >/H C i m (R 2 , C) is identified with R N for some natural number 
N. With any Helmholtz function ip defined around (xo,yo), there is associated 
a mapping 

i>: (R 2 ,(x 0) yo))- JS elm (R 2 ,C) 

defined by taking the r-jet of ip at (x, y) for each (x, y) near (xq, yo)- It is called 
the r-jet extension of ip. Moreover, to any family ^(x, y, A) : R 2 x R e — » C of 
Helmholtz functions, there corresponds a mapping 

(R 2 x R*,(x ,yo,A )) - JS elm (R 2 ,C) 

by taking the r-jet of *S>(x,y, A') at (a;, y) for each (x,y) near (xo,t/o) and pa- 
rameter A' near Ao- 

By a similar proof to that of ordinary transversality theorem/ 114|). we have: 

Lemma 6.1. Suppose a finite number of submanifolds W\, W%, . . . of Helmholtz 
r-jet space J^j elm (R 2 ,C) are given. Then, any Helmholtz function ip(x,y) de- 
fined around (xo,yo) is approximated (in C°° topology) by a Helmholtz func- 
tion ip(x,y) defined around {xo,yo) such that the r-jet extension of ip(x,y) is 
transversal to any Wi- Moreover, any Helmholtz wave ^(x, y, t) is approximated 
by a Helmholtz wave ^{x,y,t) such that j r *& is transversal to any Wi. 



3 

Helm 



(R^,C) 

ei + ie 2 , 



By using Lemma l6.ll we show Theorem l3.3l and Theorem l3.4l In Jj 

with coordinates xq, yo', a = a\ + 102; b = b\ + if>2, c = ci + ic2, e 
/ = f\ + i/2, h = h\ + U12 and k = k\ + ik 2 , we set 

^ = {0 = 0, b = 0, c = 0}, 

which is of codimension 6. Moreover, we define a submanifold W$ of codimen- 
sion 4 by the equations 



ei ci 


+ 


h 


fl 




fl Cl 




e2 c 2 




b 2 


h 




h c 2 





61 



bl 
62 



Cl 



ai 
a 2 



ei 
e 2 



0, 
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and 



bi a\ + e\ 
b 2 a 2 + e 2 _ q 

bi c 2 

together with a = 0, b\c 2 — b 2 c\ = 0, minus a locus VF 2 5 of more degenerate 
singularities, which is of codimension > 5. The definition of is from the 
idea of the iterated Jacobian Further, we set 

Wi = {a = 0, b lC2 - b 2Cl = 0} \ (W x 6 U Wl U Wg 4 ), 

which is of codimension 3, and 

Wi = {a = 0} \ (Wf U W| U VF 3 4 U 

which is of codimension 2. 

From Lemma 16.11 any Helmholtz function \p(x,y) is approximated to a 
Helmholtz function whose r-jet extension is transversal to the above subman- 
ifolds in J^ clm (R 2 ,C). This implies the following from the Whitney theory 
f |14|1. The transversality of j 3 ijj at (xo, j/o) to W§ implies that t/j has the regu- 
lar phase singularity at (xq, i/q). Similarly, the transversality to Wf implies the 
fold singularity, and the cusp singularity as a mapping from a plane to a 
plane. When a function has fold singularity, generically, there are two possibil- 
ities of phase singularities: hyperbolic and elliptic singularities. However, from 
Proposition 13.21 there is no elliptic phase singularity for Helmholtz functions. 
This shows Theorem 13. 31 Furthermore, Lemma Vr> . 1 1 claims that the transversal- 
ity theorem holds even for Helmholtz waves. Therefore, Theorem 13.41 is proved 
in the same way as above. 

Remark 6.2. The transversality to W% and means that does not 
intersect to W 2 and Wf, for the two parameter family of Helmholtz functions. 

7 Bifurcation problem of phase singularities of 
non- linear Schrodinger waves. 

We can apply our method to study phase singularities appearing in non-linear 
waves. Actually we treat local analytic waves or formal waves. Note that we 
would have to find other methods for the study of global structure of phase 
singularities of non-linear waves due to the existence of soliton solutions (see for 
instance [T7|). 
Let 

i*t + 2**x + f( x > *) = 

be a Schrodinger equation for a complex valued function ^ = ty(t, x). Here / is 
a real analytic function on R x C, for instance, f(x, = , 5| 2, I'. 



ei ci 


+ 


h 


fi 




fl C! 




e 2 c 2 




h 


h 




h c 2 
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Set \& = u + iw. Then, in the case f(x, = l^l 2 ^ the equation reads 

u xx = 2w t - 2(u 2 + w 2 )u 
w xx = -2u t - 2(u 2 + w 2 )w 

Let us denote by 

J£(R 2 ,C) :={j r *(to,x ) |i*t + ^*xx + /(af,*) = 0} 

the Schrodinger jet space. Then we see J£(R 2 , C) is a submanifold of J r (R 2 , C). 

In Jg(R 2 , C), the condition \P = gives a smooth submanifold of J£(R 2 , C) 
of codimension 2. Therefore, the phase singularities appear on the (i,a;)-plane 
at isolated points and, on the x-line, the phase singularities appear momentarily. 
The fold singularities form a submanifold of codimension 3 in J|(R 2 , C). The 
fold locus is tangent to the x-line in codimension 4. 

The condition "J" = ty x = also gives a smooth submanifold of Jg(R 2 , C) of 
codimension 4. Therefore, degenerate phase singularities appear in a generic two 
parameter family of solutions, where IP = ty x = 0. Consider the condition \P = 

= ^xx = 0. If * (t , a; ) = y x (t ,u ) = 0, then the condition ^ xx (t , x Q ) = 
is equivalent to that i^ t (to,Xo) + f(x a ,0) = 0. If /(x ,0) = 0, for instance if 
/(a;,*) = I*! 2 *, then ^ xx (t , x ) = if and only if * t (i ,a;o) = 0. Thus any 
bifurcation on t of phase singularity with ^(t ,x ) = ^ x (t , x ) = ^^(t^, x a ) = 
is degenerate (^ t {to, x o) = 0). 
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